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ON THE LOCAL UNIQUENESS PROBLEM FOR PERIODIC 
SURFACE WAVES OF PERMANENT TYPE IN A 
CHANNEL OF INFINITE DEPTH! 


D.H. Hyers and J.A. Ferling 


Introduction 


We shall be interested in studying wave motions at the surface 
of a heavy” liquid in a channel bounded by parallel vertical walls. 
we assume that the motion is the same in each vertical plane parallel 
to the walls, so that the problems considered are two dimensional, 
and that the flow is incompressible, non-viscous, and irrotational. 

There are two characteristic difficulties of such problems. First, 
the wave profile (surface) is an unknown of the problem so that even 
the region of flow is unknown in advance. Secondly, one of the con- 
ditions to be fulfilled along this unknown or “ free” boundary is 
non-linear in character. It is obtained from Bernoulli’s law in- 
volving the square of the velocity, the particle ordinate and the 
pressure, by stipulating that the pressure is constant along the 
free boundary (E.g., at an air-water interface the pressure would 
be atmospheric pressure). 

In most discussions of such problems the standard procedure is 
to linearize the equations by considering only smal] motions and 
neglecting second order quantities. One also replaces the (unknown) 
free boundary by a fixed straight line (equilibrium surface of liquid) 
in stipulating the linearized boundary condition. The resulting 
linear problem for a known region constitutes only a first order 
approximation to the original problem, but is comparatively easy 
to solve. 

One of the simplest types of such wave motion is a traveling wave 
of “permanent type.’ This denotes a disturbance which travels down 


l. Fee paper represents results obtained under contract Nonr 228-09 
NR 041/152) with the Office of Naval Research. 


2. The adjective heavy signifies that we do not neglect the force of 
gravity. 
3. Seee.g., Lamb's Hydrodynanics, chapter 9. 
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the channel at a constant velocity of propagation c, and without 
change of form. In more mathematical terms, this means that the 
motion is steady when viewed from a coordinate system attached to 
a point of the profile moving with the constant velocity c. 

In the early part of the century there was considerable discussion 
as to whether these waves of permanent type actually exist, and for 
a time Lord Rayleigh himself had doubts about the matter. A number 
of different approximate methods of solution had been devised, but 
the first mathematical demonstrations of the existence of exact solu- 
tions were given in the early 1920’s by Levi-Civita [1] and by 
Nekrasov* for the case of periodic waves in an infinitely deep channel. 

Let A denote the wave length, g the acceleration of gravity. The 
characteristic dimensionless parameter may be taken to be Y = 


——> where c is the velocity of propagation as explained above. 


Qne2 

Levi-Civita’s procedure was first to transform the problem as ex- 
plained below in section one, making use of the periodicity as well 
as of the other hypotheses to introduce a new independent (complex) 
variable { which varies over the unit circle in the (-plane. This 
has the advantage of eliminating the first characteristic difficulty 
of an unknown domain for the independent variable. At the same time 
he introduced a new dependent variable w = 7 + iG related to the 
velocity. The condition of constant pressure along the free boundary 


was transformed into the non-linear amma = -Ye~?7 sin 6 to 


be fulfilled along the circumference { = e'? of the unit circle in 
the (-plane. By assuming an expansion for both & and 6 in terms of 
a certain real parameter 44, he was able to determine the coefficients 
in both expansions and finally to demonstrate the convergence of 
these expansions in the neighborhood of “ = 0. Thus the existence 
of exact solutions was demonstrated locally. 

However, it seems that no adequate discussion of uniqueness of 
solutions of this problem has yet appeared in the literature. This 
paper represents a first step in this direction, being concerned 
only with local uniqueness. 

We use a method explained by 0.P.Friedrichs [2] to reformulate 
the problem as a problem of “bifurcation.” (cf. sec. 2.) In this 
approach, the parameter of Levi-Civita arises naturally, and Levi- 
Civita’s ad hoc assumption of the expansions in powers of yu is 


avoided. 
As a by-product of our method, we prove that a conjecture of 


Levi-Civita®, namely that every solution with Y > 1 is simply related 
to a solution with 5 < 1, is true at least “‘locally”, i.e. for Y in 
the neighborhood of each integer. Our results are given in section 3. 


4. We have not had access to Nekrasov’s papers [5], but only to an 
stract of them given by Levi-Civita. 


5. ai. p. 284, 
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1. Levi-Civita’s formulation of the problem. 
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Figure 1. 


Let c denote the velocity of propagation of the wave and let ? 
denote the wave-length(assumed independent of the time). The (x,y)- 
coordinate axes are to be 1t.ixed to the profile | of the wave, witt 
the origin at a point of l and with the y-axis pointing vertically 
downwards. Our assumptions imply that the flow when viewed fron 
this coordinate system is steady, two dimensional, irrotational, and 
incompressible, and is periodic with the period A. 

We denote the components of the particle velocity relative to the 
coordinate axes by (u,v). The actual velocity (relative to an 
observer fixed to the earth) is then (u-c, v). We shall assume 
that the actual velocity is zero at infinite depth, so that u-+e¢ 
and v — 0 as y ~ wo Setting w =u - iv, we see that w = u + it 
is the (relative) velocity vector in the complex z-plane, where 


+ ty. 
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The irrotationality of the motion and the incompressibility of 
the fluid allow us to .:..roduce the velocity potential ?, the stream 


function ¥, and the complex potential X = + where 
z 


The functions ® and ¥ being determined only up to additive con- 
stants, we may adjust the latter so that ®= ¥ = 0 at the origin. 
Thus in particular ¥ = 0 on the free boundary /| and x(0) = 0. 

The region of flow (including the free boundary |) in the z-plane 
maps into the half-plane #2 0 in the x-plane (Figure 1). 


From the assumption that the function w(z) has the real period x 
it follows easily that x (z+) = x(z) +2. Thus if we look upon 
the complex velocity w as a function of x, w will have the real 
period cA. For this reason it is convenient to introduce a new in- 
independent variable ¢ defined by 


Considering now the mapping from the z-plane to the C-plane we see 
that it js many to one. However two values of z differing by A cor- 
respond to the same point © under this mapping, so that w is a single 
valued function of . Clearly, points of the free boundary t cor- 
respond to points of the circle | ¢| = 1 and the “‘point’’ y = + © cor- 
responds to ¢ = 0. If we restrict ourselves to a period strip in 
the z-plane, bounded by the profile and the equipotential surfaces 
= ‘2c, we have a simple map onto the unit disc in the [-plane 
(see Figure 1). 

The origin in the z-plane is mapped into the point © = 1, but so 
far z = 0 is an arbitrary point of the profile 1. For the present, we 
leave this choice of origin abritrary, but fixed. Later on, the 
choice of z = 0 will be made® by requiring that dw /dt shall be real 
and negative for ¢ = 0, which has the effect of specifying one di- 
drection in the C-plane. 

To obtain the boundary condition along the free surface | we 
observe that the pressure there is constant. Thenby Bernoulli’s law, 


(1.1) Y(u2 + v2) - gy = const. along l, 


where g denotes the acceleration of gravity. 
In order to simplify the problem we introduce a new dependent 
variable w = t + 16 in accordance with the equation 


where the principal value of the logarithm is used. We shall always 
assume, with Levi-Civita, that for all points of the flow 


6. See [4] p. 101. 


w = Log w/c, 
— | 
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(1.3) 


c 


Hence the function w({) given by (1.2) is analytic fo 
In terms of the new variable w = t + 16, the c 
becomes 


(1.4) 


We set 5 = pe*?, so that vy = 2n O/cX is a dimension! 
potential. Differentiating (1.4) with respect to 7 and simplifying, 
we get the boundary condition 


dt 
(1.5) — = -Ye~** gin 0 for |¢| = 1 
dep 


where y is the dimensionless parameter: 
(1.5) Y= re / 


The original problem of finding a steady flow w(z) with period 4 
satisfying the boundary condition (1.1) along 1 and the condition 
w = c at infinity is now replaced by the problem of finding, for a 
given positive number Y, a function w(C) satisfying the following 
conditions: 


(a.) = t+ 10 is analytic for |¢| <1; 
(b.) wld) is continuous for |¢| < 1; 
(c.) w(0) = 0; 


(d.) = sin for |¢| = 1. 


For convenience, we shall refer to this problem as Levi-Civita’s 


problem. 
When the profile 1 of the wave is symmetric about some vertical 


line the wave is called symmetric. If the y-axis is chosen to be 
along this vertical line, it can be shown that the condition 


(1.8) = w(Z) 


must hold for a symmetric wave. Conversely if (1.8) holds, it is 
obvious that the wave is symmetric with the y-axis of as an axis 
of symmetry. 


7. For further details concerning the relationship of this problem to 
the original problem see [1], or tal, pp. 97-101. 
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2. Another formulation of the problem. 


Following K. O. Friedrichs [2], p. 160-178, the above problem will 
be formulated in a somewhat different way. However, we do not assume,as 
Friedrichs does, that the solutions are symmetric in the sense that 
= wl 

In the classical linearization of Airy the non-linear boundary condi- 

dt 

tion = - Ye 6 is approximated by the nonear = 

and this leads to the approximate solution w(%) = 40" where Y is the 


positive integer n. 
In trying to find true solutions ‘‘near’’ these approximate solutions, 
it is convenient to rewrite the principal boundary condition (1.5) in 


the form 
d 

(2. 1) + nd = nb - sin 6 Citi = 2}, 


Our next objective is to formulate the problem in terms of an inte- 
gral equation. To this end we study first the simpler problem 


dT 
(2. 2) — bly) along |f| = 1, 
fo) 


where for the present, b (~) is a given function. Later we shall replace 
b(p) by the right hand side of equation (2.1). 

Let H, denote the space of all functions af) which are analytic for 
[Z| < 1, continuous for |%| £1, and vanish at { = 0. With real numbers 
as scalars, H, is a vector space, and with the norm: 


H, is a Banach space, as is easily verified. 
Let H, denote the subclass of H) for which the two quantities 


d . 
lim 2. Re [wloe'?)] and — Re [wle'?)] 
p=1 dep 


exist and are equal for -7 £ g< 7. 

We then let B denote the class of all real valued functions b(¢) on 
-7 £ @ <7 which are continuous and satisfy the two conditions b(77) = 
b(-77) and f* b(p)dp = 0. With the usual maximum norm, B is a Banach 


space 
LEMMA 2.1 For each b € B, the function 


(2.3) 

is analytic for |{| < 1 and ImF(pe **) tends to bla) as p tends to one 
(o <1). 


67 


ON THE LOCAL UNIQUENESS PROBLEM 


Proof: If we put [=e t@ and perform some easy calculations we find that 
(1 - p*)b(o) 
-7 1 + p* - 2pcos(e-d) 


TmF(C) = % ble) dy += J” 


The lemma now follows from the known properties of Poisson’s integral, 
since the first integral is zero by hypothesis. 

Let w(Z) belong to H, and satisfy the boundary condition (2.2), where 
C= pe*? and b(¢) is - element of B. Then observe that 


dT 


Now consider the function 


(2.5) G(t) =F(U) + (07%) 


where F ({) is defined by (2.3). By lemma 2.1, G'Z) is analytic for lZ|<1 
and its imaginary part has a limit as P ~ 1 - which vanishes identically 


on the circle |%| = 1. Hence G(Z%) is constant for \t| < 1, and since 
G(0) = 0 the value of the constant is zero. Therefore 

d d 


and Fi) =~ . Hence Fg) = 0. Now by (2.3) the Taylor 
expansion of F(Z) about [= 0 is 

-7 


Thus the condition F ("0 ) = 0 is equivalent to 


To obtain an integral representation for w(Z), we divide equation 
(2.6) by ("*? and integrate, obtaining 


(2.9) wat at* 


Then we substitute (2.7) into (2.9), making use of (2.8). The integra- 
tion with respect to Z may be carried out term by term and thie result is 
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tmp 


7 


n-1 imp 


m=1 


blple*"? log (1- dp + 


where the principal value of the logarithm is used. From (2.10) it 
follows that p = (0)/n! 

These considerations establish the following result. 
LEMMA 2.2: Let 6 be an element of B. Then if w() is an element of Hy, 
which satisfies (2.2), it may be written in the form (2.10) or (2.11), 
where 44 is a constant, and b must satisfy (2.8). 


In the following section we shall also need the converse of lemma 2.2 
for the case n= 1. For n= 1 equation (2.11) becomes 


LEMMA 2.3: In the case n = 1, if 6 € B satisfies (2.8) and alZ) is de- 
fined by (2.12), then w belongs to H, and satisfies the boundary con- 
dition (2.2) 


This lemma may be verified by straight forward but somewhat tedious 
calculations (see [2], pp. 188-178). 


3. Theorems of existence and uniqueness. 


In this section we study the non-linear integral equation obtained 
from (2.10) by replacing 6() by the right member of (2.1), and obtain 
a locally unique solution depending on the two parameters Y and 4. Next, 
the cases where n> 1 are reduced to the basic case n= 1. For this 
case, 4 may be taken as real without loss of generality. After estab- 
lishing three properties of this two parameter family or solutions, it 
is shown that the auxiliary condition arising from (2.8) is satisfied if 
and only if a certain relationship holds between ¥ and 4. This leads us 
back to the solution of Levi-Cicita’s problem. 


LEMMA 3.1: The transformation b ~M [b] defined by 


(3.1) (Le?) bi 


n 
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K (z) + = z"log(1- 2), 

is a linear bounded transformation of the Banach space B into the Banach 
space I1,, for each positive integer n. 


Proof: Clearly Vv [b] lies in H, and is linear in 5 for b €B. To prove 


the boundedness, observe that 


cm /2 


0 


log |sinx| dx = - ~ log 2 


and that the principal value of log(1- e# (%9)) is 


5 \a-p) + log |sin 


Setting 6 = ped in (3.1) and using the above results we have 


log 2 + sgn(¢p - a)} + 


+ (2 + 4m) log 2 + 3m}, 


THEOREM 3.1: Let ‘3 
Ye sin Glo), 


) + i6(~) =wle'?) and put ht) = " where n is a positive 


where T\o 
integer. To each n there corresponds a positive number 0) such that, for 
every real ¥ with y-nl < o, and every complex there is at most one 
solution of the equation 


9) M [b (w,a)] + ub 
noon 


n 
lying in the & neighborhood of the origin of the space H.. If tn ad- 


dition (4) < on 2, then there is exactly one solution w = w(t.) iL) in 


this neighborhood. 


Proof; The partial derivatives of the function b, = n6~-Ye-%7sin 6 with 
respect to 7 and 6 can be made arbitrarily small by taking 7, ©, and 
|¥-n\ sufficiently small. Also by lemma 3.1, the transformation b — 
M [b! is linear and bounded. Hence there exists a positive number 0, 


such that 


whenever the numbers ly- nl, are each less than The state- 


ment of uniqueness follows easily from this Lipschitz condition. 
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To prove the existence, we first note that equation (3.2) is 
satisfied by taking = 0, =n, and w = 0. Also, when < 0/2, 
l|uh || < 6/2 and it follows by using the Lipschitz condition 
with {2 = 0 that the transformation w — M[b (w,¥)] + wh, takes the 
neighborhood || || < 0, into itself. Hence we may use the method of suc- 
cessive approximations to establish the existence of the solution in 
this neighborhood. 8 

The next theorem shows that we may restrict ourselves to the case nel, 


THEOREM 3.2: Using the notations of the preceding theorem, we have 


for each real a in the neighborhood of Y = n and each complex « in the 


neighborhood of zero. 


Proof: For a given value of n and k < 1, let 5 = 0, be chosen in ac- 
cordance with theorem 3.1, and let € = min (0;,0,). Now choose ¥ so that 
< € and choose so that < €/2, Putting = /n, we have 
~ 1| < €/n < 83, and | <3 (1- kk). Hence by theorem 3.1 there is a 
unique solution w = @;(C, ¥, 4“) of equation (3.2) for the case n= ], 


Writing out this equation (cf. (3.1)) and replacing 2 by %" we get 


a \ \ 
rn(j+1) 
dp+uC, 


J 


where w, = 7,(p) +24, (p) has period 27. The integrand being per- 

lodic, we may replace the limits by -n7 and n7 if we simultaneously 

divide by n. Then we are ready to make the change of variable ~ = ng, 

and obtain: 

nj 


nb, [w (ng ,)n,¥,] dy, + ut. 


jz=i 


A change of the index of summation by putting m = nj +n gives 


nb, [w, (nq, ut®, 
in which the index m of summation varies over the proper multiples 2n, 
3n. ... of n. However, since b, b, lw, 1s continuous and 
has the period 27, it follows that 6,(nq@ is orthogonal to the functions 
e™'? in case m is not a multiple of n. Thus these other values of m 
may also be included in the summation, and it follows from (3.1) that 
¥/n, u) is a solution of (3.2) in which 


8. See [3], pp. 133-134. 
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b n@, (ng) - Ye (np) sin G, (np) nb, [w,(np),6 


4 
But this solution is unique and the theorem follows. 


From now on we limit ourselves to the basic case n=], so that Y will 
be selected in the neighborhood of Y = 1. For this case = w (0) and 
by a proper choice of the origin in the z-plane, we may take 4 real with- 
out any loss of generality.? 


THEOREM 3.3: The solution w= all,¥,u) of the equation 
(3. 3) M,[b,(w,y)] 


whose existence was established in theorem (3.1) has the following pro- 
perties. For every real value of ‘Y and of yu in suitable neighborhoods 
of y=landu-0, 


(i) w(l,¥,) has derivatives of all orders with respect to 5 andy; 


(ii) ¥,u) where @ begongs to H, and has deriva- 
tives of all orders with respect to a and 4; 


(iti) ¥, = wll, uw), so that all thesolutions are symmetric. 


Proof: Since M, [6] is bounded and linear, and since b, lw, Y) = 0-Ye~ 77 
sin 6 is indefinitely differentiable in 6, 7, a, it is easily seen that 
the right member of (3.3), looked upon as an element of Ho, has Frechet 
differentials of all orders with respect to the triple (@,¥, 4). More- 
over, each of these differentials is continuous in all of its variables 
uniformly with respect to the increments $a, dy for 
lsy| < 1, Su] <1. Consequently, it can be shown that the right mem- 
ber of (3.3) is of class C (n , in the sense of Hildebrandt and Graves 
[3], for each positive integer n, and the property (i) now follows by 
applying their theorem IV, p. 150, 19 

To establish properties (ii) and (iii) we define a sequence of ap- 
proximations by 


(3. 4) = 0; = ye om sin 6] + uh, 


where + Trivially, = 0 for = 0. Making the induction 
assumption that = 0 for = 0, it follows from (3.4) that = 0 
for = 0. Hence the limit 


wll, ¥, w) = lima lt, 


taken in the sense of the norm in Ho, also vanishes for u=0. From this 
result together with property (i), property (ii) is easily established. 


9. See [4], p. 101. 


10. One can also establish property (i) by using the method of successive approx- 
imations directly. 
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Property (iii) can be proved inductively in a similar manner, if we 
observe thafw,(¢) = @ (C) implies that O, is odd and 7,is even in 9, 
so that -37 
b,(p) =Q - Ye ~ sin & 


is odd. 

Our next objective is to show that the auxiliary condition (2.8) for 
n= 1 can be satisfied if and only if a certain relationship holds be- 
tween the parameters and Y. Substituting the solution y, of 
(3.3) into (2.8), where b = 6 - Ye~47sin 6, gives the relation !1 


(3.5) Viy,u)= sing - ye-37'9) sin 6 (p)} dp=0. 


In order to show that this equation can be solved for Y near pu = 0, 
we first divide by 4, obtaining the equation 


A 


(3.6) =f" sing {6 - yen 9, 


Now the function w(Z, y, KL) satisfies the equation 


from which it follows that @(Z, 1, 0) = h (0) = ¢. Differentiating 
V(y, u) under the integral sign with respect to y and setting y = 1 and 
0, we get 
” sin - sin’y dy =-7 

Clearly equation (3.6) is satisfied for a = 1 and » = 0. Hence by the 
ordinary implicit function theorem, (3.6) may be solved uniquely for y 
in terms of 4, in the form Y = ¥\(y), this function being defined near 
ft = 0. We summarize our results in the following theorem. 


THEOREM 3.4: There exist positive numbers 4, n such that for each value 
of the real parameter 1 = w'(0) satisfying lu| <n, there exists a uni- 
, solution y, w of Levi-Civita’s problem satisfying |1 - y| < 6, 
llw|| < 8 This solution is symmetric in the sense that w{Z) = a(f) 


Proof: To establish the uniqueness, let y, @ be a solution of Levi- 
Civita’s problem. Then since w = 0 for ¢ = 0, it follows by the mean 


value theorem for harmonic functions that i Op) dp = 0. Consequently 


1l. Since has the symmetry property (iii), the function b - 
is odd in Y, so that relation (2.8) reduces to (3.5). 


| 
| 
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b € B, where b = 6 - ye~%7 sin 6. By lemma 2.2, w@ satisfies (2. 12) 
or equivalently (3.3) and 6 satisfies (2.8). By theorem 3.1 there is a 
unique solution of (3.3) providing y and 4 are in the neighborhood of 
Y = 1 and uw = O and & is in the neighborhood of the origin in the space 
H). Since b satisfies (2.8), y and yu satisfy (3.6), so that y= ylu) is 
uniquely determined by 4. 

The existence is proved by reversing this argument, beginning with 
the solution @ = @4(¢, y, 4) given by theorem 3.1, and using lemma 2.3 
in place of lemma 2.2. 

Let the choice of z = 0 on the profile h be made so that 4 is real 
and negative, as described in section one. Since & is symmetric by the- 
orem 3.3, it follows that the origin z = 0 is at a crest or else at a 
trough of the wave. We next show that the origin z = 0 is located at a 
crest when uu < 0. 

Putting 7, = Re at 1), we see that since the solution w/{) is deter- 
mined by the parameter 4, 7) is a function of u alone in the neighbor- 
of « = 0. Now using theorem 3.3 with y = y{u), we see that T) 1s a con- 
tinuously differentiable function of alone, and that T | ), 
where 7 <4) = Re &(1) is also continuously differentiable in 4. It fol- 
lows that T, (0) = F, (0) = 1, since w = { for u = 0. Since 7(0) = 0, we 
see that 7) 1s negative since 4 is negative. Thus w = ce” is real and 
less than c at z = 0, which therefore cannot be a trough and must be a 
crest. 12 

The argument just concluded in the last paragraph also shows that we 
may use 7 or the velocity w. = ce’° at a crest as the basic parameter 
in place of yu, tf desired, in stating a uniqueness theorem analogous to 
theorem 3.4, 


THEOREM 3.5: There exist positive numbers 0, € such that for each & 
satisfying 0 <1- y < 8, there exists a unique solution w = ail) of 
Levi-Civita’s problem satisfying the inequalities < € and w'(0) <0, 


Proof: By theorem 3.1 for the case n = 1, there exists a 5 > 0 such 
that equation (3.3) has a unique solution w = w,(%, a, uu) satisfying 


ly- 1) <8, lul <% 8, lloll < 


and by theorem 3.3, w,(%, y, 4) has derivatives of all orders with re- 
spect to % and u. By expanding wilt, y, 4) in powers of Y and yu by Tay- 
lor’s formula with remainder, determining the first few coefficients in 
this expansion by substitution into (3.3) and then substituting the re- 
sulting expansion into the auxiliary equation (2.8), Friedrichs ([2], 
pp. 179-182) has shown that (2.8) or equivalently (3.6) reduces to }3 


12. See Oj, section 9. 
13, In equation (4), p. 182 of [2], the coefficient of c? should be -1. Equa- 
tion (3.8) also agrees with Levi-Civita’s relations (18), p. 312 of [1]. 
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(3.8) 


It follows that (3.8) can be solved uniquely for 4? in terms of y, and 
since by hypotheses u = w'(0) < 0, for u in terms of y. More precisely, 
there exist positive numbers 5’ < 5, and € < %5, such that (3.8)has a 
unique solution = under the restrictions 0 < 1-y< 6’ and 0 < 
- uw < €. By lemma 2.3 the function w= @,(0, y, uly)) is a solution of 
Levi-Civita’s problem. By restricting Y even more if necessary, say 
0<1-~y< 8 < 8’, we can insure that lee Since lw" (0)| < lew it 
follows that - uly) = -w'(0) < €. Hence our solution is unique under 
the conditions of the theorem. 

In terms of our original quantities, putting U - u-c, } 
have: 
THEOREM 3.6: There exist positive numbers 6 and € such that, for every 
value of the wavelength \ and the propagation velocity c satisfying | < 
1 - (gh/2rc?) < 8, there exists a unique non-trivial wave of permanent 
type in an infinitely deep channel whose (actual) particle vélocity vec- 
tor (U,V) satisfies U*? + V2 < €2. There are no such waves satisfying 
this restriction on the velocity and 


/-~v, we 


0 < (gi - 1< 8, 


except the trivial one for which U =V=0. 


REFERENCES 


1, T. Levi-Civita. Determination rigoureuse des ondes permanentes d’ampleur 
finie. Mathematische Annalen, V. 93 (1925) pp. 964-344. 


K.O. Friedrichs. Functional Analysis and Applications. Lecture Notes by 
F.A. Ficken, New York University, 1949-50. 


. T.H. Hildebrandt and L.M. Graves. Implicit Functions and their Differentials 
in General Analysis. Transactions of the American Mathematical Society. 
V. 29 (1927) pp. 127-153. 


. J.B. Serrin and G.E. Latta. Notes on Hydrodynamics. Part I. Princeton Uni- 
versity, 1951-52. 


5. A.J. Nekrasov. On Stationary Waves. Bulletin of the Polytechnic Institute of 
Ivanovo-Vosnesenskii, V. 3 (1921) pp. 52-65 and V. 6 (1922) pp. 155-171 (In 


Russian). 


University of Southern California. 


A NOTE CONCERNING HOMOGENEOUS POLYNOMIALS 


J. C. Morelock and N. C. Perry 


1. Introduction 


In classical algebraic geometry there has been considerable 
interest in specific instances of homogeneous polynomial surfaces 
F(z,, Xo, Xa, ay? = 0 of prime degree p, which are invariant under 
the cyclic transformation 7 defined by: 


For example, it is a known result, obtained apparently by actual 
trial, that the homogeneous polynomial of degree 5 with the largest 
number of terms which is invariant under JT is 


This paper proposes to give a general treatment of these instances 
by characterizing al! such polynomials for any prime p > 3, by a 
method which employs elementary number theory only. 

We note, to beginwith, that a surface F(x,, x,, x,) = 6, Se 
invariant if each of its terms x{ xt x$ we is transformed by T into 
the expression E* x4 where s is of the same residue class 
(modulus p) for all terms. This situation will be referred to here- 


after as S-invariance. 


2. Necessary and sufficient conditions for s-invartance: 


Required conditions on a, b, c, andd are: 
b +2c + 3d = s'+ kp (Os s' ¢p-1) 


Since a, 6, c, and d are non-negative integers where each is 
at most equal to p, it is clear on inspecting the last equation 
that k can take on only the values 0, 1, 2, and 3. 

Diophantine solutions of the two equations are easily obtained 
by elementary methods as: 


* L. Godeaux’s Geometrie Algebrigue 
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a=n 
b=r-2n-s' + (2-k)p 
c=n-2r+s' + (R-1)p 
If k = 0 we have the four inequalities: 


(1) OSn<p 


(2) OF r-Qn-s’ +2p 
(3) 0 
(4) 


1A 


n-2r+s‘' - p 


O<crcp 


From inequality (3) we have 2r - s’ + p <n and from (2) it follows 
r-s' +2 
that n < eee, Thus n must lie within the range defined by 


2pt+r-s' 
Since from (3) 4r -2s° +9p 2nwe 
may replace 2n in (2) by 4r - 2s’ + 2p obtaining 0 < r - 4r + 2s’ 


2r-s' 


- 2p - s‘ + 2p which simplifies into 3r < s‘. 

Hence, to obtain an invariant term of degree p, necessary re- 
restrictions on r and n are (for the case k=0)_ first that r be 
chosen on the range 0 < r < s'/3, and then n be chosen on the range 
2r -s' r 

The argument by which these restrictions on" and r were obtained 
from inequalities (2) and (3) is clearly reversible. Hence, since 
the obtained solutions to the Diophantine equations are the only 
possible solutions, the conditional inequalities 0 <r < s'‘/3, 
2r-s' +psns 4(2p+ r- s'), constitute with the equations 
determining 5 and c, a set of necessary and suficient conditions 
for an invariant term of degree p, for the case k = 0. 

By similar arguments, necessary and sufficient conditions can be 
obtained for the cases, k = 1, k = 2, and k = 3. 

We summarize all these results in the following: 


Theorem I: Necessary and sufficient conditions that a polynomial 
term x¢ «5 xf + of degree p = 3 shall be s-invar- 
tant under T are - 

loa=n, b=r=-2-s' + (2 —k)p, 
c=n-2r+s' + (R-1)p, d=r 
(where k=0, 1, 2, or 3) 
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2. If k=0, O<rcs'/3, 2r-s' 
If k=1, O¢rds'/2, 


4(s' +1) <r < 1/3 (p+ s'), 2 -s' S nS -s) 


If k=2, OS n SK(r - 


YU(p+s' +1)<r<1/3 (Qp+s'), 2 ) 


If k=3,r =p, n=0 =s' = 9) 


3. The number of terms in the most general invariant 
polynomial for any s’. 


By an intricate analysis the possible choices of the pair (n,r) 
satisfying Theorem I, can be enumerated for any s’. The method employed 
[5], uses ordinary number theory only, but involves an elaborate 
bread-down into cases where all possible combinations are considered 
of (r-s) even or odd, p of the form 6a +1, s’ of the form 6/, 
68+1,68+2, 56 +3, for each of the first three sets of condi- 
tional inequalities in Theorem Il. 

We note that for a given s’ distince pairs (n,r) produce different 
polynomial terms, because n is the exponent of x, and r is the ex- 
ponent of x,. Then since any two of the inequality sets are non-over- 
lappong in the least one of the variables mn andr, any duplication of 
polynomial terms in the enumeration is impossible. The same pair 
(n,r) may arise for different s’ but in this instance the exponents 
of x, and Xs will be different, since 6 and c are strictly monotonic 
functions of s’. 

The results obtained can be summarized - 


Theorem II: For each s‘ such that 1< s' <p-1_ there are 


p? +6p +11 : 
P invariant polynomial terms. For s' = 0, 
there its one additional term as noted under the case 
+ 6p + 17 
k = 3, or terms. 


6 


These results check with the classical theorem stating the number 
of terms in the entire homogeneous polynomial* of degree p, for: 


(p? + 6p +11) (+ + 2)@ + 3) 
6 6 


+ 6p +17 
6 


+ (p -1) 


* Snyder and Sisam, Analytic Geometry of Space. 
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The method described in this paper permits a classification of 
the general homogeneous polynomial into p sets, where to each set 
there corresponds an invariant surface. Any subset of any one of 
these sets also yields an invariant surface. 

Thus, the approach used makes it possible to form invariant poly- 
nomials without trial and error, and also affords anew way of counting 
the terms inthe most general homogeneous polynomial of prime degree p. 


4. Conclusion; A numerical example. 


We illustrate the application of Theorem I and II for the case 
p = 5. Geometric usage of such developments can be found in Dessart 
[1], and Hutcherson [3] and [4], for polynomials of degree 5, 7, and 11. 

For p = 5 the inequalities of Theorem I yield the following 
tabulation of pairs (n,r): 


| «2 


or 
@ 


Ww 


oOo 
oF 
o % 


7 


oo 


r 
0 
0 

1 
0 
1 
2 
3 
5 


As described in Theorem II this tabulation yields a 12-term 
polynomial for the case s’ = 0, and 1]-term polynomials for each 
of the cases s’ = 1,2,3, and 4. Thus, in all, we have 12 + 44 or 
56 terms inthe most general homogeneous polynomial, as described by 
the general expression (p + 1)(p + 2)(p + 3) /6 when evaluated for 

5. 

These five polynomials are listed below as computed through use 

of the (n,r) tabulation and the expressions for band c in Theorem IT: 


§ 3 2 2 ie 4 3 5 
+ t+ + Xy + Xy 
4 4 2.2 2 3 2 


ae 2 2 2 we 3 4 
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k n k 
5 0 | 
0, 4,2 
0 
0 
2 0,1 2 
1 
3 0 
Fy; 
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4 23 4 3 


2 2.2 2 2 2.9 4 
4 3 4 2 
3 3 4 
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Certain Non-factorable Polynomials 
Doyne Holder 


Teachers of algebra in secondary schools find it convenient to 
be able to determine quickly if a trinomial of the general type 
ax* + bx +c is or is not factorable. A proof is give: here that 
a certain type is non-factorable. 

Theorem: In the general trinomial ax? + bx +c, if a, 6, ande are 
all odd integers, then there are no linear integral factors 
of ax? + br +c. 


Proof: If ax* + bx +c is factorable, let these factors be such that 


+ bx +c = (mx +r)(nx +s). 
Then 

(1) 

(2) re=e 

(3) «ms +nr= 


Since both a and ¢ are odd, it is trivia! that m, n, r, and s are 
all odd, hence all non-zero. Furthermore, the products, ms and nr, 
are both odd - each being the product of two odd numbers. 

By (3) ms +nr= 6, and 6 is odd. But the sum of two odd 
numbers cannot be odd; hence in the case where a, 6, and c are all 
odd there are no linear integral factors of ax? 4 bx 4c. 


The above theorem can be extended to polynomials of higher degree 
by altering the proof only slightly as shown here. 


Theorem: If = +a. + tax ta, po lynomia | 


whose coefficients are integers and if a,, c,, anda (1) 
are ali odd, then there are no linear integral factors of 
a(x). 
Proof: If there is a linear integral factor of a(x) let it be 
d(x) =m +k 
It is a trivial that m must divide a, and k must divide a,. Since 
both a, and a, are odd, then both m and k are odd, hence both non- 


zero. 
Furthermore d(1) divides a{l) and as a({1) is odd, d(1) must be odd. 
But d(l)=m +k , botim and k being odd, and the sum of two 
odd numbers cannot be odd. 
Therefore, if a,, a, and a(l) are all odd, there can be no 
linear integral factors of a(x). 


o 
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THE LUCKY NUMBER THEOREM 


D. Hawkins and W. E. Briggs 


The lucky numbers of Ulam resemble prime numbers in their apparent 
distribution among the natural numbers and with respect to the kind of 
sieve that generates them [1]. It is therefore of interest to investi- 
gate their properties, bearing in mind the analogies to prime number 
theory. 

The lucky numbers are defined by the following sieve. If S, is an in- 
fmite sequence of natural numbers t, , (m = 1,2,3, ...), one obtains 
for n > 1, from S, by removing every for which divides 
m. S, is the sequence 2,3,5,7,9, ... of the number 2 followed by the 
odd integers in increasing order of magnitude. S, is the sequence of 
natural numbers. The sequence of lucky numbers is S = Lim S|, that is, 


2,3,7,9, 13, 15,21, ... . This definition differs trivially from Ulam’s. 


_ Two properties are basic for the investigation of saymptotic proper- 
ties of lucky numbers. By the definition if s, represents the m-th lucky 
number, then 


(1) =t,, foralla<s. 


This follows from the fact that ae is the first number that will 
be removed from S, in forming S,,,, and is at the same time less than 
any number to be removed later on. Also by the definition, if R(n,x) is 
the number of numbers not greater than x inS,, and [x] denotes the 
greatest integer in x, then 


_R(n - 1, x) 


s 


(2) R(n,x) =R(n - 1, x) n 22, 


n-1 


This fundamental recurrence relation has the following solution, in 
which {x} denotes the fractional part of x and 


o, = (le 1/2)(1 - 3)(1- ... (1- Vs,_,), 


R(i - 1, x) 


(3) R(n,x) = Le] o + 


n 2 2, 
im2 
= [x] o + E(n,x). 
(4) 0 (nx) en. 
The properties is S now develop by a series of stages. The first is 
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to find bounds for s,. If one puts R(n,s)..) =n +r, which by (1) one 


may do for 0 <r then 


n+r 


n 
(5) 


by letting r =0. From (5) 


Putting p, = Yo. 


Summing from 2 to n - 1, one obtains ~, ~ p, > L which implies 2 > 
log n, or t=2 


(6) o 


n ~ log n’ 


In =sne¢r, (0 Sre n), one may set r=n andr =n — l, 


since s, 2 2n for n > 2, which gives, by (3) 
2n = E(n,s,) 
Hence, by (4) and (6), for n > 2, 


(7) s, >nlogn Son-1 > (2 1) logan. 


2n 


It is now possible to show that the remainder term E(n,x) is o(n) when 
x = s,. For from (1) and (3), 


(8) R(n, + E(n,s ) 


Let a(n) be the integer defined by Sa(n) 2 n and Sa(n)-1 <™ By (7) 


(9) a(n) < 3n > n> No» 


log n 
Now split the sum E(n,s,) into two parts E, and E,, where 


(10) R= 1, 5) 


| 
n+1i l o 
| 
> n¢+1 1 
n- 1 
| 
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o RG-1, s) 
(11) 

Clearly 
(12) 
log n 
In E,, because of (1), put all R(i - 1, s,) = Rlaln), s,) =n, so that 
by (7) 
n 
n nlog log n 


(13) 
i log i log n 


i=3M/logn 
It is now possible to write for n 2 2, 


oO 
n 


(14) 


Again using the substitution P,= Yo, , one obtains 


1 
"P, n 


By summation /, = logn + 0 (log n) and, therefore 


log n 


and from this, (1), and (3) 


n + o(n) 


On 

(15) is the analogue of the Merten’s theorem for prime numbers, and 
(16) is the analogue of the prime number theorem. 

These results, especially (16), confirm a conjecture of one of the 
authors based on stochastic arguments (2). They support the observation 
that the asymptotic distribution of prime numbers is not, except in de- 
tails, a consequence of their primality, but characteristic of a wide 
class of sieve-generated sequences, of which the lucky numbers are an 


example. 
By using the results recursively in (14), S. Chowla has shown that 


the asymptotic value of (15) can be improved to 


(17) s, =n logn +5 (log log n)* + o[n(log log n)?). 


Since the corresponding result for prime numbers (where P,, is the n-th 
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prime number) is 


(18) p= n logn +n log log n + o(n log log n), 


it follows, with only a finite number of exceptions, that s, >p,. With 
necessary calculations, this presumably will confirm Ulam’s conjecture 
Ss, > P, for all n. 
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CRITICAL CURVES IN SEISMIC EXPLORATION 


H. G. Helfenstein 


1. Foreword 


In seismic prospecting a charge is exploded at a point O on the 
surface of the earth and the time for the seismic waves to return 
to the surface is measured in various points in the neighborhood 
of O. The velocity of sof the waves in the uppermost layer down to the 
reflecting discontinuity surface D can be measured separately. By 
means of these data we can plot points of the “‘time-distance\surface” 
defined as follows. 

Introduce a system of coordinates with O as origin, Z-axis 
pointing vertically downwards, any x,y plane coinciding with the 
surface df the earth. At every point F = (u,v,0) we measure a value 
of the travel time T(u,v) and a total distance w = sT travelled by 
the wave. The surface with equation w = sT=!)(u,v) is the “time - 
distance surface’’S. 

It is easy to see that the fictitious surface S obtained from 
our measurements and the physical surface D are connected by a con- 
tact transformation T [1]. According to the law of reflection the 
ellipsoid E of revolution with foci 0 = (0,0,0) and F= (u,v,o) and 
major axis w #!l(u,v) is tangent to D. The equation of E is given by 


(1) w =Vx2 4 y24 224+ x)? + (v y)24 22, 


where u,v,w =fQ(u,v) are kept constant. If on the other hand we 
keep x,y,z fixed and vary u,v,w equation (1) represents one sheet 
of a hyperboloid H of revolution. The surface D is the envelopeof 
all the surfaces E obtained by letting u,v,w vary on S, and vice 
versa S is the envelope of all surfacesH obtained by letting x,y,z 
vary on D. 

We have shown [2] how this remark can be used to determine graph- 
ically points of D together with their tangent planes from a finite 
number of points of S. Our present problem is concerned with the 
position of the recording instruments. We shall show that the data 
gained from certain critical curves may be insufficient for the de- 
termination of the transforms of these curves under 7. In practice 
the instruments will be distributed over a whole area around 0, and 
hence these critical curves are mostly of no practical importance. 
It is an amusing geometrical problem, however, to give a complete 
characterisation of them; they turn out to be arcs on certain 
hyperbolas. 
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2. Transformation of Curves 


We first turn to the following question: What information about 
D can we obtain if we measure the reflection times along a certain 
curve? In orther words, assuming that we know a curve K on S which 
surfaces are then possible discontinuity surfaces? 

Three close points of K are transformed by T into three close 
ellipsoids E, and their common point (or points) tend to a limiting 
position which will be a point of a curve k, image of K under T. All 
possible surfaces J) must contain k. 

We find k analytically according to the theory of envelopes by 
substituting the three functions u(t),v(t),w(t) representing K in 
(1) and differentiating twice with respect to t. It is more convenient 


first to rewrite (1) in terms of the quantity r = Var + ye 22, 


This yields the equation 
(2) ux + vy — wr - Yu? + v2? + w*) = 0 


It is linear in the new coordinates x,y,r, and slightly more general 
than (1). The latter equation defines a real surface E only for 
u2 + v2 < w2, while (2) represents for u2 + v2 > w2 hyperboloids or 
planes. 

We call the separating cone u“ + v? = w* the critical cone. Our 
surface S is situated entirely in its interior. From (1) we also 
conclude that the tangent planes of H (which coincide with those of 
S) are never inclined by more than 45° towards the u,v plane. Hence 
no curve on S makes an angle greater than 45° with the u,v plane. 

Diferentiating now (2) with respect to t we obtain the following 
linear system for x(t),y(t),r(t) defining the curve k. 


u(t) x + v(t) y- w(t) r =%(4? v? - 


2 


2 


a“ 


3. Definition of Critical Curves. 


The solutions os system (3) remain undetermined if its augmented 
matrix 


Yiu2+y2- 
(A) %(u2 + v2 w2)’ 
Y(u2 + v2 


has a rank less than 3. A curve for which this condition holds in 
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every point will be called critical. 

The determination of the image of a critical curve under T requires 
higher derivatives than the second. It is therefore not advisable 
to place the seismometers on the projection ontotheu,v plane of a 
critical curve because we might by chance pick up by our measurements 
just such a critical curve, 


4, Determinatton of all critical Curves. 


The four determinants of order 3 in matrix (4) are all Wronskian 
determinants. Their identical vanishing entails the existence ofa 
constant matrix (A;,) each of whose rows contains a non-zero vector 
such that the following four equations hold identically in t. 


(5) = 0 
(7) + ut Az ,(w? - C*) = 
(8) Ag + (w? - (7) = 
where =Y v2 , 


We observe that every straight line passing through the origin 
and every straight line making an angle of 45° with the u,v plane 
satisfy these conditions. Hence they are critical. 

Since each of the equations (5-8) represents either a plane or a 
surface of the second degree passing through O all the other critical 
curves are conics passing through O. Let us first assume that the 
plane of such a conic is not perpendicular to the u,v-plane (all 
the results remain valid in the latter case however). This plane 
is then given by an equation 


(9) w= Aju + 


and the projection of the conic on the u,v-plane by 


(10) au? + buy + cv? + du + ev = 0. 


Substitution of (9) in (6) gives a relation between u and v which 
must be the same as (10). Hence a constant A # 0 must exist such 
that 


MA,,=4, AAy >= 
2 
6A4 -1) =a, = 6b, 


87 


88 MATHEMATICS MAGAZINE (Sept.-Oct. 


Since we ruled out the straight lines at least one of the numbers 
a, 6, c does not vanish. Therefore, Ang # 0. Hence the last three 
equations yield 


e2fi-D, b=2frr, ¢ = fl -D, 


where f is a non-vanishing constant. The same solutions are obtained 
by intersection of (9) and (7, 8). It is easy to see that these 
necessary conditions are also sufficient, and our critica! conics 
can therefore be represented by 


dy, - Iu’ + 2A ,A,uv (2 - 1)v? + du + Ev =0, 


where A,, A,, D, E are arbitrary parameters. 

Some elementary analytical geometry leads from this expression 
to the following complete geometrical characterization of the cri- 
tical curves: The set of critical curves consists of 

a) the straight lines passing through the origin, 

b) the straight lines making an angle of 45° with the u,v-plane, 


c) the ellipses passing through 0, whose minor axes are parallel 
to the u,v-plane and between whose eccentricity € and whose 
angle w of their plane with the u,v-plane, the following re- 
lation holds: 

(11) e = V2 sina, 


1) the parabolas ; 
make angles of 45° with the u,v-plane, 


e) the hyperbolas passing through O whose transverse axes are 
parallel to the u,v-plane and for which (11) also holds. 


5. Projection of the critical curves on the u,v-plane. 
For the projections of these critical curves we find the following 
properties: 


Every conic C’ of the u,v-plane passing through O is the projection 
of a critical curve. Its plane is found by the following rules: Its 
intersection with the u,v-plane passes through O and is 


a) parallel] to the minor axis of C’ if C’ is an ellipse; 
b) parallel to the transverse axis of C’ if C’ is a hyperbola; 
c) perpendicular to the axis of C’ if C’ is a parabola. 


In all these three cases its angle with the u,v-plane is found from 


(12) e'= tang 
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where €'is the eccentricity of C’. 

Not the whole of these conics, however, is ‘‘dangerous”, there 
are only certain arcs which should be avoided. It is clear that 
those arcs of a critical curve which are in the exterior of the 
critical cone w? = u? + v?, and those arcs making angles greater 
than 45° with the u,v-plane, can never be obtained by actua! measure- 
ments, since their images on D are imaginary. 

“Forbidden” are therefore only those arcs of a conic C’ for which 
the corresponding arcs of C are in the interior of the critical cone 
and which at the same time make angles less than 45° with the u,v-plane. 

According to (11) the planes of critical ellipses and parabolas 
are either outside or tangent to the critical cone. Hence no are 
of an ellipse or parabola is forbidden. 

The plane of a hyperbola C however intersects the critical cone 
in two real generating lines g,, g, which make an angle / given by 


(13) tan = Vv-cos 2w . 


From (11) the angle be seen the asymptotes of C can be computed, and 
comparing with (13) it is found that / = a. Using (12) and the fact 
that the transverse axis of C is parallel to the u,v-plane, it can 
easily be proved that this relation is preserved in the projection. 
Thus 6’ = a', 28" being the angle between the projections By &, of 
By» 6 and 2a' being the angle between the asymptotes of C'. Hence, 
if the two perpendicular lines g;, g4 are drawn from O to the asy- 
mptotes of C’, they define four angular sectors; the arcs Ct of C’ 
included in those two sectors whose bissector is perpendicular to 
the transverse axis of C’ are projections of arcs of C situated 
within the critical cone. 

Let us determine now those arcs C, of C whose tangents make angles 
with the u,v-plane less than 45°. 

Since the transverse axis of C is parallel to the urv-plane there 
are four of these arcs, situated symmetrically with respect to both 
axes and extending to infinity. We determine their boundary points 


in the following way: 


If a line-element is situated in a plane making an angle w with 
the u,v-plane and if its projection includes an angle x with the 
trace of its plane then its slope is determined by 


tang T = sin X tang @, 
From <45° we conclude, using (12), that 
sin X < l/e’ = cos a’, 


where a’ is the acute angle between an asymptote and the transverse 
axis of C'. Hence the limit points of the dangerous arcs of C’ are 
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characterized by 


(14) 


where Xo is the angle of the tangent of C’ with its transverse axis. 
This equation expresses the fact that the tangents in the boundary 
points are perpendicular to the asymptotes, and therefore parallel 
to the above mentioned lines ae 85° 

Since the slope of a tangent of a hyperbola is numerically greater 
than the slope of the asymptotes we have always 


(15) >a’. 


Hence if a’ > AS” there is no value X o which satisfies sim |taneously 
(14) and (15) and therefore there are no forbidden ares on the hy- 
perbola C’. 

If a’ < 45 , however, it is always possible to draw 4 tangents 
t; to C’ perpendicular to the asymptotes. Their points of contact 
A, B, C, D are the boundary points of those 4 arcs C, of C’ which 
we are looking for. 

Summing up, we may state: 


a) all straight lines in the u,v-plane passing through O are 
forbidden; 

b) every hyperbola passing through O whose asymptotes include 
an angle 2a’ less than 90° (measured in that sector which 
contains the curve) contains one or two forbidden arcs. They 
are the common parts of the arcs C’, and C', determined above. 


Given an arbitrary hyperbola C‘ with asymptotes making an angle 
less than 90° we can convice ourselves of the truth of the above 
statement by letting O vary along C’. The points A, B, C, D are 
independent of the position of 0. Let E be the intersection of the 
parallel line to t, through B with C’. If C' is given by € = a cosh t, 
y= b sinh t, the parameter values of A and E are found from tanh t,= 
- b2/a? and te = 3t,. Our figure illustrates the case when O is 
situated between A and E. Other cases to be considered are: O be- 
tween A and B, O = A, and O outside of BE. ze 
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University of Alberta 


Dear Mr. James: 


“The paper presented is the work of my twelve year old boy who 
took more than a year to gather courage to submit it for editorial scrutiny. 

You may be interested to know that when my son first received a sub- 
scription to Mathematics Magazine about two and a half years ago, he was 
aghast to note that he couldn’t understand a single thing in it. With each 
successive issue, however, his understanding unfolded (he is a self-taught 
mathematician) until now he awaits each issue with eagerness, and recent- 
ly was able to submit his solution to one of the Proposals published in 
your last issue. He considers his subscription to Mathematics Magazine 
one of the finest presents he ever received!’’ 


Sincerely yours, 


Sylvia Bergman 


91 
, 
q 
+hea 
if, 
y 


N 
> 
2 
> 
~ 
8 
3 
— 
3 
= 
> 


See 
LAPD JEDoms 
| 


SOLVING DIFFERENTIAL EQUATIONS 
WITHOUT COMPLEX NUMBERS 
Michael J. Pascual 


The only place in a course in real differential equations that 
complex numbers arise (embarrassingly, I’m afraid) is in the study 
of the linear equation with constant coefficients. This seems to 
imply that the theory is not complete in itself, without trespassing 
into the complex number field. The purpose of this paper is to show 
that one need not introduce complex quantities at all to obtain the 
solution even when the auxiliary equation does not have real roots. 

The equation we are concerned with is 


where p and q are real constants. For convenience this can be put 
into the form 
(1) y" -2ay’ + (a? + b*)y = 0 


where a= - and b =” |q - (p2/4)|. We first solve (1) when a =0 


dv 
by setting y'=v, y" = ae that (1) becomes 


dv 2 
(2) v — + b*y =0 
dy 


which is easily solved and yields solutions 


t by? 


the double sign inside corresponding to the one in (2). The solutions 
for y are 


y =A sin (B+ bx) and Inlby +’ B+ b*y*) =+ bx +A 


where A and B are arbitrary; the first solution corresponding to 
the minus sign of (2) and the second solution corresponding to the 
plus sign. With a little manipulation these can be put into the form 


(3) y = cy cos bx + cy sin bx 


(4) y= +c 


in which the ambigious sign is taken care of by the arbitrary character 
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of c, and c,. 
Now we are ready for the general case with a # 0 and b # 0. The 
change of variable y = ue** transforms (1) into 


(5) u” + b4u = 0, 
which is the same as (2): hence 


u =c,cos bx + C > sinbx 


so that 
(6) y = e°*(c,cos bx + c,sin bx) 


(7) y = + ce” 


where (6) holds for the plus sign in (5) and (7) holds for the minus 
sign. 
Finally for the only case we have not covered, namely when b = 0 
we solve (5) easily to obtain 


u= + Co 
so that 
(8) y = + cp) 


Now that all possibilities have been covered in a perfectly natural 
manner without introducing complex quantities, one might for expediency 
show that by solving the auxiliary equation 


m? ~2am + (a2 + 67)=0 


and “examining’’ the roots we can immediately set down the solutions 
(3), (4), (6), (7), or (8) depending upon whether the roots are pure 
imaginary, real and negatives of one another, complex conjugates, 
real and unequal in absolute value, or real and equal respectively. Of 
course once (1) has been disposed of then, bythe usual development 
of operators in factored form, higher ordered equationscould be 


solved. 


Siena College, New York. 


or 
u =c,e* 


SOLUTION OF TRIANGLES ON THE SLIDE RULE 


Vv. 


C. Harris 


The following examples show how the law of sines may be used 
with the slide rule to solve triangles in case three sides, or two 
sides and the included angle, are given. Methods of solution by use 
of the law of sines in other cases are well known. Hence, subject 
to the limitations of the slide rule, the law of sines suffices to 
solve any triangle. 

The method involves successive approximations. The adjustment 
made at each step is estimated on the basis of the difference between 
180° and the sum of the angles obtained by the previous approxi- 
mation. Skill can be developed rapidly and the triangle usually 
solved in three or four settings. 

We take up first the case of three sides given. 

Example 1. Given a = 4.8; b = 7.3 and c = 8.2. Opposite the 

largest side (8.2) on the D scale set 99° on the S 
scale. (We will call the scales D and S hereafter and these will be 
the only scales with which we shall be involved.) The angles on S 
cpposite the three given sides on D and their total are: 


Cc: 90° 
B + 63° 


Total 188.9° 


Since the total is not 180°, the angles obtained do not yield a 
solution. The triangle is not a right triangle since if it were, the 
angles just found would have to be the solution. Moreover the tri- 
angle must be acute since the sum is more than 180° and if C were 
increased, both A and B would be decreased but (as the slide rule 
shows) by an amount less than the increase in C and the total of 
the angles would always be too large. We now try to reduce the sum 
of the angles to 180°. If we should take 8.9° off angle C so that 
opposite 8.2 on D we would have 81.1° on S, then we would also be 
taking off of angles A and B and our result would look like this: 


Cc: 
B: 61.7° 
A: 35.3° 


178.1° 


Total 
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We would have compensated too much and wish instead to allow for 
this. We should take off less than the discrepancy and we estimate 
that we should take off about 8° instead of 8.9°. The result is thus: 


A: $5.4 


Total 179.2° 


We have still taken off too much, but we see that changing C by 8° 
took off only 1.2° fromB and .5° fromA. Hence most of the .8° 
which the total of the angles needs to make 180° will come from 
increasing C . If we try 82.6° we find 


82.6° 
B : 
A: 35.5° 


Total 180.0° 


Hence these values give the solution. 


The work may be put down compactly as follows: 
First Second Third 
setting setting setting 


82 82.6 
61.8 61.9 
35.4 35.5 


179.2 180.0 


In case the largest side does not start with the largest digit, 
it is convenient to multiply or divide the sides by the same integer 
to have the largest side have the largest initial digit. The purpose 
of this step is to make all readings occur to the left of the index 
of S so that the slide need not be reversed. Also, if the sum of 
the angles is less than 180° when the first setting is made, then 
the largest angle of the triangle is obtuse, and the supplement of 
the angle on S must be used in evaluating that largest angle. 
Example 2. Given a = 131; b = 77 and c = 92. Divide the sides by 

2 and solve the similar triangle as indicated: 
First Second Third Fourth 
setting setting setting setting 


100 101 101.4 
35.4 35.2 35.1 
43.8 43.6 43.5 


179.2 179.8 180.0 


a 


Cx 90 
63 
‘ a-=4.8 A: 35.9 
188.9 
65.5 A: 90 
B: 36 
45 C: 44.6 
170.6 
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The angles obtained in the fourth setting are the solution. When the 
slide is moved to the right, angles B and C are diminished in the 
present example. Hence A must be increased to obtain a total of 180°. 
The largest angle must be obtuse and is the supplement of the angle 
read on S. Note that in increasing A, since angles B and C decrease, 
we must increase A by more than the discrepancy in each trial. For 
example, the discrepancy at the first setting in 9.4°,so0 we es- 
timated an adjustment of 10° which we then applied, but this turned 
out to be too small. 

When two sides and the included angle are given, the procedure 
is much the same. Again, if the larger side does not start with the 
larger digit, it is convenient to solve a similar triangle first 
and then multiply or divide by the proper number to get the unknown 
side. 


Example 3. Given 6 = .91; ¢ = 1.25 and A = 71°. First, multiply 
the sides by 5 and solve the similar triangle 6 = 4.55; 


c = 4.25 and A = 71°. The work follows: 


First Secon Third Fourth 
setting setting setting setting 


71 71 7] 7] 
46.8 40.5 42.4 42 
90 43 68 67 


207.8 174.5 181.4 180.0 


Side a is read on D opposite 71° on S as 6.42, so in the origina 
triangle we have a=1.28. Note that we still let the largest un- 
known angle be 90° for a first approximation. In using 63° for C in 
the second approximation we show that not very much skill was being 
employed since almost no allowance was made for the expected decrease 
in B. Consequently the total change in the sum of the angles turned 
out to be 5.5° too much. 


In a few instances, some angles must be read on the ST scale. 
Also, the number of settings needed varies. However, the preceding 
remarks and a little experience will show the user what happens to 
the small angle or angles as the largest one is adjusted. This en- 
ables one to vary the adjustment of the large angle in the correct 
direction and cuts down on the number of settings needed. 

When compared to alternative methods, the above method seems 
simple in idea and rapid in practice. It also has the advantage of 
using the same trigonometric law as that commonly employed in solving 
triangles under the other cases, a law which can be explained to 
students with little mathematical background. 


San Diego State College 


a= A: 
b= 4.55 -B: 


A NUMBER SYSTEM WITH AN 
IRRATIONAL BASE 


George Bergman 


The reader is probably familiar with the binary system and the deci- 
mal system and probably understands the basis for any others of that 
type, such as the trinary or duodecimal. However, I have developed a 
system that is based, not on an integer, or even a rational number, byt 
on the irrational number r (tau), otherwise known as the ‘‘golden sec- 
tion”, approximately 1.618033989 in value, and equal to (1+ V5) /2. 

In order to understand this system, one must comprehend two peculi- 
arities of the number r. They are based on tau’s distincive property 1 
that 


(a.) Take any approximation (A,) of r. Taking the reciprocal, we get a 
number (a) that is proportionately the were distance from 1/ras A 1 was 
from, but arithmetically nearer. Adding 1,? we get a number (A>) that 
is atieabuhemabaly nearer 7 than a; was to l/r but arithmetically just 
as near. Since a; is arithmetically nearer than A,, A,is nearer in both 
respects to r than A,. Repeating the process of taking the reciprocal 
and adding 1, we approach +. Now, taking 1 as A,, and expressing our 
approximations of 1 (i.e. AyA>,A3, etc.) as fractions, we get 

1235813 


Taking either the numerators or the denominators, we get what is 
know as the Fibonacci Series, each term of which is formed by adding 
the two previous terms 3. for 


1 

1 n 
( 


80 


(We designate the n“ term of the Fibonacci Series by f,» setting 3 1, 
f, = 1, This practice shall be used throughout the aticie.) 


(b.) Any integral power of r can be expressed in the form r” = AT+ B, 
where AandB are integers and, in fact, numbers in the Fibonacci Series. 
The explanation of this startling fact is really rather simple: 
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Since 


r+ =1r +0 and r2 = lr +1, 


and since A 
gt, (ie +0) + (lr +1) 4 


In the same way 
4 = (Ir + 1) + (27 + + 


In general 


~ 
~ 

+ 


4 
tht (see Note 2.) 


Can this be applied to negative powers of +? We don’t know any Fibon- 
acci numbers before 1, but it is easy to see how we can find them: 
Taking 1 and 1 as our first two, we can see that the term before must 
be 0, since that is the only number which, when added to 1 gives 1. In 
the same way, the number before that must be 1, since lis the only num- 
ber that, when added to 0 gives 1; and the next term must be ~-1, since 
no other number gives (0 when added to 1. Continuing this process, we 
get 0, 1, -1, 2, -3, 5, -8, 13, -21... . Obviously, this is alternately 
+1 and —1 times the corresponding Fibonacci numbers. But can this be 
proved to be true in all cases? It can by induction. The rule we want 
to prove, expressed as an equation, is: 


fy = 


Let us assume it true for y = 1, 2, ... n. Now by the basic property of 
the Fibonacci Series: 


f _(n41) = (-1) it f,) 
f_(n+1) = (-1)"*? fast 


The inductive proof is completed by the examples already cited. 


5 5 


Applying this to powers of r, we make a list of them from r~°tor”: 


NUMBER SYSTEM 


Now, at last, we shall get back to our concept of a system based onr. 
Like the binary system, it can have only two symbols: 1 and 0. But, un- 
like the binary system, it has the rule (2): 100 = 0115 (place the 
decimal point anywhere - it’s a general rule). But how do we find the 
numbers? We know that 1 is r° or 1.0. Next, looking at the table of 
powers of fr, one notices that 


+0 and +2 
Adding them together, one gets 
ri472 = +0) + + 2) = 2, 


Therefore, 2 <10.01 (in this system). Of course, because of rule (2), 
this can also be expressed as 1.11, © 10.0011, 10.001011, 1.10101], 

etc., but 10.01 is what I call the simplest form (that form in which 
there are no two l’s in succession, and which, therefore, cannot be 
acted upon by the reverse of rule (2), called simplification (1l= 100). 


To convert a number to its simplest form, repeatedly simplify the left- 
most pair of consecutive l’s. 

To continue with our “translation” of numbers into this system, we 
notice (after a careful examination of the table) that r* = lr + 1 and 
r~* = -lr +2, and adding them together r?4772 =3, and so 3 is 
100.01 in this system. What about 4? Well, since r 
+17 °(101.01) must equal 4, since r 1, Can this method of adding be 
used for other numbers? The answer is “yes”; just convert the number 
into the form in which there is a zero in the units column and place a 
lin it. If the method of conversion is not obvious, use this method: 


a. Change to the simplest form. 


b. If there is no 1 in the units column, you are finished. If there 
is, look in the r~? column (there can’t be any in the r~! column be- 
cause it is in its simplest form and there is a 1 in the colum next 
to it); if there is = 0 there, expand? the 1 into the r~! and r~2 
columns (that’s all); if there is 1, look in the r~* column; if there 
is a zero there,expand the 1 in the r ~2 column into the r~ and r-4 
column and the ] in the units column into the r~! and r~2 columns. If 
there is a l, look in the r~® colum; if there is a zero, expand the 
1 in the r~ colum into the r~° and r~° columns, the 1 in the 77? 


100 
Sr -8 Or +1 
re +50 = + 0 
Se 43 
| 
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column into the r~%and r~*column, and the 1 in the units colum 
into the and +r ~* columns; if on the other hand there is a 
look in the r~®column, etc. If it is the endless fraction 1.0 1010101..., 
change it to 10.000000 ... . 
We can now construct a table of integers in this system. Here are 
those from 0 to 14 (in their simplest forms) 


0-0 5 — 1000. 1001 10 — 10100, 0101 
1-1 6 — 1010.0001 11 — 10101.0101 
2- 10.01 7 — 10000.0001 12 — 100000. 101001 
3 - 100.01 8 — 10001.0001 13 — 100010.001001 
4~101.01 9-—10010.0101 14—100100.110110 
Examples of the basic processes 

1. Change 100101.111001 (equals 16) to the simplest form. The first 
pair (farthest left) is in the units andr~! column, so we simplify it 
into the r+ column, giving 100110.011001. This time the pair farthest 
to the left is the one we have just created with our new 1 in ther? 
column (the result of our simplification), added to the 1] already in 
the r * column. This we simplify into a 1 in the r* column, which gives 
us 101000.011001. Finally, we change the last remaining pair (in the 
r-2 and r-3 colums) into a 1 in the 7-1 colum, arriving at our final 


answer: 101000. 100001 


2. Change 101.01 (4) to a form with a zero in the units column. (i.e. 
a form to which 1 can be added). One can see that it is already in its 
simplest form. However, there is a 1 in the units colum, and we must 
remove it. The first thing we do is look in the r ~? column; since there 
is a l there, we look in the r~ colum. This is empty, and so we ex- 
pand the 1] in the r ~* column into the r “3 andr~* columns, getting 
101.0011. Now that the r~* colum is empty, we can expand the unit into 
@ pair in the r-2 and 7-1 columns, getting 100.1111, which has a zero 
in the units colum. We can now add 1] to it: 


100, 1111+1=101. 1111 =110. 0111 =1000. 0111 =1000. 1001 (five) 


The Arithmetic Operations 


The arithmetical operations, although they are basically the same as 
in any other system, are, in practice, quite different because of the 
peculiarities of this system. As our first step in all of them, we eli- 
minate zeros, which would only hinder us, and show the place values of 
l’s by actural placement in columns. For instance, four (101.01) would 
be [1] |] [1], the heavy line representing the “decimal point”. The 
necessity for this step results from the fact that, though in the sys- 
tems to which we are accustomed, the steps of addition are simple enough 
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to be performed mentally, this is not so in the tau system; nor is each 
column non-dependent on the one to the left of it. It is thus neces- 
sary to have lined columns in which to carry out the work. 

Now for the actual processes, we shall start with addition. The example 


10010.0101 / 9\ would be represented " If a | 
+ 1010.0001 \+6 1 1 


In this set-up it can be seen that we have a pair, consisting of a l in 
the r> column and one in the r* column. This we simplify into a 1 in the 


r> column. 
1 1 | 


Now, however, we have no obvious way to continue. We are left with two 
l’s_ in the same column. We can neither add them together to give 2 (as 
we would in the decimal system), nor is there any simple “carrying” 
operation. We must, therefore, change this to a form not having two l’s 
in the same column. We will start by expanding one of the 1’s in the r’ 


ila | 
ai ala tat | 


column: 


Now we can simplify the pair we have just created in the r/ and units ; 


columns: 
1/4} fala 1 
1 


A 
Aj 
We shall now use the same type procedure for the 1’s in the r~* colum. 
We expand one of the 1’s there: in |x 1\2 1 
| 
1 


and the one in the r~! and r~ columns: 


| 


and simplify in the r~* and r ~° columns: 


| if | 
A | AVA 
and express our answer in ordinary form, writing l’s in columns with an 
un-crossed-out 1 and 0’s in the columns where all have been crossed out: 


100101.001001 (15) 
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For general rules as to procedure, I believe that these will do 
in most cases, (These general rules and the ones for the other 
processes are not the types of rules that, if disobeyed, would give the 
wrong answer, but merely guides to the quickest way to get the right 
one): 


a) Expand only when that is the only way to remove a 1 from the same 
column as another, regardless of whether this will result in the same 
situation in another cclumn, but only if no more simplification can be 


done. 


b) Simplify whenever possible, and, if there are two or more pairs, 
always simplify the one farthest to the left first. Because of this rule, 
simplification should never result in two 1’s in the same colum, i.e, 


(1412) should be simplified into lata tat rather than | | 


Subtraction is the next process I shall describe. As in addition, we 
set up the numbers in columns, but here we shall assign negative values 
to the 1’s from the subtrahend. For instance, to find (11-6) we set up 


| 1 1 | } 
We now ‘‘cancel” the 1 and the -1 in the r~ column, giving 


1 1 ] 


Subtract ton 


1 
-1 


1 


-1 


Next, we expand the 1] in the r 2 column, getti 
1 


-l 
Again, we cancel, this time in the r/ column, 


1 


1 -A 
And again we expand, this time the 1 in ther’ column, getting 

1 1 A 
-1 | | 


| 
| 

ng 

-A 
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For the third time we cancel, (in the r? column) giving 
A|A | 1 1] 
-A 
This we treat just as we would if we were adding and arrived at this 
stage; by expanding one of the l1’s in the units column we get 


Next we simplify the pair in the units and 7~! columns getting 
Aj A| AJ AJ A 1 A 


1 


and the pair we thereby form in the r } and 


r? column getting 1 | 
ja 


Finally, we expand one of the 1’s in the 
colum, getting 


A ‘| 
| alala 
columns) Ge 


A 
-A 
or 1000.1001 (5) 1 A 1 


and simplify the resulting pair @n the r~ 
our final answer 


For subtraction it is harder to formulate a general rule, but I think 
it would suffice to say: Cancel whenever possible and simplify or expand 
whenever that would permit cancellation (also remember not to confuse a 
1 with a -1 ({*{1|-1|/4]1/2| -2|) and not to make the mistake of ‘‘expand- 
ing’ a -1 into two+l’s.) After all -1’s have been removed by cancella- 
tion, proceed as you would with an addition example. 


Multiplication 


Multiplication involves nothing new. We simply place the partial pro- 
ducts as we do in the decimal system, and add. For instance: 
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101.01 
x 100.01 (3 x 4) 


Setting up the partial products, we get 10101 or 
10 


101 
| 
yi y ll | 
and (from now on it is simple addition) expanding one of the l’s 


We now simplify the pair in the units and r~! column 


Big: 
1 


and the pair we thus produce: 


I 
A 
A 


] A 


and again the pair this simplification produces, getting: 


| alata | 
2 column, and one of the 


11 

1] | 

Finally, we simplify the pair in the r~? and r~? columns, and then 
the one inthe +" andr~° columns, giving our final answer: 


aa || 
1jA Ail 1 | 

or 100000.101001 
Division 


Next we expand one of the l’s in the 77 


2 


Division is quite different in this system, and is, in fact, rather 
odd. The only things it has in common with ordinary division are the 
basic principles behind it, the way the example looks, and the movement 
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of the ‘‘decimal point’’ to eliminate any figures to the right of it in the 
divisor. It is best explained ” an example: 12 divided by 2, or 


which, after i brn point’’, is 


| 


Now, since we are dividing by i | [1], if there are anywhere two 1’s 
with two spaces between them (the spaces can be empty or full), they can 
be crossed out and a 1 placed in the quotient above the rightmost of the 
two. This crossing out in no way signifies that the 1’s should not be 
there, but is merely equivalent to, in long division (decimal) the subtrac- 
tion of the product of the number placed in the quotient and the divisor 
from the dividend. Since the number placed in the quotient can only be 1 
(placed in any column, of course), we merely subtract the dividend (placed 
in that same column), i.e. cross it out, It is obvious that once the whole 
dividend has been crossed out, the group of 1’s in the quotient, after 
being changed to the simplest form, will be the complete quotient. Getting 
back to our original problem, we see that we do have just such a set of 
1’s in the r~! and r~* column, and so we cross it off and place a 1 in the 


| | 


But now, you may say, there are no more pairs of 1’s spaced in that way; 
what shall we do? The answer is our old pair of friends, expansion and 
simplification. Since they do not change the value of a number, if either 
of those processes yields a set of 1’s spaced correctly, that set can be 
crossed off and a 1 placed in the quotient just as though that set were 
part of the original number. Since in our problem it is so far impossible 
to simplify, we shall expand. iy the 1 in the r’ column, we get 


quotient, getting 


No such set yet. However, when we expand the 1 we’ve just placed in 


1 {al 


the r> column, giving 


| 
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we have not one but two such sets (remember that all we need is two 1’s 
with that certain separation, regardless of intervening and crossed out 
1’s), one made of the 1’s in the r~ and P columns, and the other of the 
1’s in the r? and r® columns. Crossing them both out and placing the 1’s 
in the correct places in the quotient, we get 


| 1 1 
/ 


and since there are no more 1’s in the dividend, our number in the quotient 
is the complete quotient, and so 1010.0001, or 6, is our answer. This 
time the general rule is: Always take that course of action that will place 
your next 1 (i.e., a 1 in the quotient—set in the dividend) farthest to the 
left. By a course of action, I mean a series of expansions and simplifi- 
cations and the exchange of a set for a 1 in the quotient that follows; or 
simply that exchange, if the set is already there. (I did not obey this rule 
in my demonstration so that I could show the process in a simpler way.) 
This is so that the answer be in its simplest form. 

By the way, the processes of addition, subtraction, and often multipli- 
cation, can be performed together by writing the addends, the subtrahends, 
and the partial products in one set of columns; for instance: 2x3+4+3-5: 


10.01 
x 100.01 + 101.01 + 100.01 -1000.1001 


and working it out: 


| 
1 
1 partial products 1 | | 
i| {a} ja 
1 1 addends | 
-1} | |  subtrahend bri, | | 
0 


Now that we know these four processes, we have a much better way of 
finding a number in this system than merely repeatedly adding 1’s until 
we reach it. For instance, to find thirty-seven, we can multiply 6x6 and 
add 1; to check the arithmetic, we multiply 7x5 and add 2: 


1010.0001 (to check) 1000.1001 
x 1010.001+1 x 10000.0001 + 10.01 
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What’s more, we are not only able to find integers, but, since we can 
divide, we ought to be able to find fractions also. Let us try. First, we 
shall attempt to find 4%. We begin by setting up our division: 


Our first set, as can be seen, wil] have its leftmost 1 in the r! column. 
We therefore expand the 1 in the r~ column: 
| 


| 


The other one we need to complete the set is a 1 in the r~ 
we get by expanding the 1 in the units column: 


wee 
After we exchange our set for a 1 in the quotient (r~2 column), we notice 
that our remainder is 1 (in the r~* column), Since 1 is the number we 
started with, the next figure in the quotient and the next remainder should 
be the same as these. However, the question is, where in the quotient 
shall we place it? Since our first 1 was in the +~ column (because we 
moved the decimal point) and our remainder is three places to the right of 
it, in the r~! column, our next 1 in the quotient should be three places 
to the right of the first 1 there. Since the next remainder will bear the 
same relationship to the first remainder as the first did to our original 1, 
the following 1 in the quotient wil! be three places to the right of our 
second 1. Since this can be carried on indefinitely, it appears that 4% 
expressed in the Tau System is .01001001001 (any ‘‘doubting 
Thomases”’ may carry it out a few places to see). 

Before we go on to other fractions, it would be wise to mention some- 
thing about 1 in the Tau System. As you can easily see, 1=.11=.1011= 
.101011=.10101011 etc. It is, therefore, equal to the endless ‘‘fraction’’ 
.10101010.... (just as in the decimal system 1=.9999999....). If we can 
now take this fraction and expand the leftmost 1, and then expand the 
1 in the r73 column, so as to prevent the occurrence of two 1’s in the 
same column, and then expand the 1 in the r~” column so that there are 
not two 1’s in that column, etc., we will get .01111111.... If, on the other 
hand, we start by expanding 1’s in other columns, we get: .100111111...., 


column; this 
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-10100111111....., 101010011111....., etc. Therefore, if you multiply 
-01001001001.... (4%) by 10.01 (2) and get .1001111111...., this does not 
mean that 2x,=a fraction, but merely shows a different way of rep- 


resenting 1. 
To get back to fractions, we can make a list of them just as we did of 


integers before: 
1/2 = .010010010010 
1/3 = .001010000010100000101000 
1/4 = .001000001000001000001000 
1/5 - .00010010101001001000001001010100100100 


1/10 = .000010000100010100001010001010101000100101000001001000100 


00000001... 
Of course, finding these fractions is immeasurably harder than finding 
%, and with 1/10 I had to work it out 5 or 10 times before I got the correct 
answer, as there is much room for error. 

By the way, no fraction can be terminating in this system, since that 
would mean that it could be expressed as the sum of a group of integral 
powers of Tau. Since all the powers of Tau can be expressed as the sum 
of an integer and an integral multiple of Tau; if the integral multiples 
‘“‘cancel’’ (e.g. 7-3-3745 = 1) the result will be 
an integer, and if they don’t (e.g. 27-3-—3r+5=r+2), it will naturally be 
irrational. However, when we have an endless series, this paradox is 
detoured by admitting the fact that Lim AT+B with A and B always 
integral can be a rational fraction if A and B > ~, 

The Tau System has a good many other interesting and unusual char- 
acteristics, and investigation by the readers of some, such as the frequency, 
occurrence, and nature of numbers with a 1 in the units column (when in 
simplest form) might prove interesting. I do not know of any useful application 
for systems such as this, except as a mental exercise and pastime, though 
it may be of some service in algebraic number theory. For instance, the 
numbers expressible in the Tau System in terminating form consist of 
all the algebraic integers in R(/5), and some of the properties of numbers 
in this and other systems might correspond to facts about associated fields. 


Definitions Invented for Work in the Tau System 


Expand: alter three successive figures of a number by changing ...100... 
to ...011... The result is the same in value as the original, because 
of rule (1). This does no¢ mean change zeros to ones and ones to 
zeros; just this specific change. 

a ‘fy: the reverse of expand; alter the figures thus: change ee O11... 


.-100... One speaks of simplifying the i’s in the r"~* and r” 
ye into the r”™ column. Also, one speaks of expanding the 1 in 
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the r” column into the r™l and r"-2 columns. 


Simplest form: that form of a number which has been simplified until no 
more simplification is possible. It therefore has no two 1’s in suc- 
cession. It also has the fewest 1’s and is the easiest form to work with. 

Columns: just as in our decimal system we speak of a units column, a 
ten’s column, a hundred’s column, a tenth’s column, etc., in the Tau 
system we speak of a units column, ar column, a r~! column, ete. 

Pair: two 1’s in succession. 

Cancellation: a change of the form 


-1| 
Set: in division, two or more 1’s arranged with the same spacing as the 


1’s in the divisor (regardless of intervening 1’s). A set can be “‘ex- 
changed’’ for a 1 in the quotient. 


1a1so true of -1/r; there are other numbers which have similar properties, 
e.g. there is a number S between 1 and 2 for which S?=S?+S+1. Ed. 


2because 


3this is the basic property defining the Fibonacci Series. 

4There is also a more complex proof which involves multiplying the ex- 
pressions like 2r+1 by r, giving 2r7+r, and expanding. r* into r+1. 

5This is a restatement of r?=r7—147"2. 

6By changing the 1 in the r column to 1.1. 


Tf you come across words (like ‘‘expand’’) used in an unfamiliar way, 
look for them in the list of definitions at the end of this article. I have 
put there all words which I have had to invent or alter for use in this 
system, so as not to break up the text by explaining them. 
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PROBLEMS AND QUESTIONS 
Edited by 


Robert E. Horton, Los Angeles City College 


Readers of this department are invited to submit for solution 
problems believed to be new and subject matter questions that may arise 
in study, in research, or in extra-academic situations. Proposals should 
be accompanied by solutions, when available, and by any information 
that will assist the editor. Ordinarily, problems inwell-known textbooks 
should not be submitted. 

Solutions should be submitted on separate, signed sheets. Figures 
should be drawn in India ink and twice the size desired for reproduction. 

Send all communications for this department to Robert E. Horton, 
Los Angeles City College, 855 North Vermont Ave., Los Angeles 929, 
California. 


PROPOSALS 


320. Proposed by M. S. Krick, Albright College, Pennsylvania. 


Find three consecutive odd integers less than 10,000 which are 
divisible by 27, 25, and 7, respectively. 


321. Proposed by Rk. 1. Coffman, Kichland, Washington. 


Given a line of unit length and an acute angle ©, construct with 
compasses and straight-edge a plane figure having an area equal to 


6 
J cos26 dé. 
0 
322. Proposed by D. A. Steinberg, Livermore, California. 


Show that if x is any real number, a any non-zero real number and p 
any non-negative integer, then — 


j k 
111 
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323. Proposed by Leo Moser, University of Alberta, Canada. 


Some problems which mention no integer explicitly have unique integer 
solutions. For example, «What is the smallest perfect number?* has the 
solution, 6. Prove that for every positive integer n, there is a problem 
which mentions no integer explicitly, but which has n as its unique 
solution. 

324. Proposed by J. M. Howell, Los Angeles City College. 

Prove that; 

cos x/(ltsin x) + log (ltx) =1 + O(x4), O<x<1. 
325. Proposed by R. B. Kiltie, Montclair, N. J. 
Given the family of concentric ellipses 
x2/A2 + A2y2/M2 = 


(1) Show that each has the area 7M. 
(2) Find the points at which the ellipses are tangent to the 


equilateral hyperbolas 4x2y? = M?. 


(3) Show that through every point for which 4x2y? < h? there pass 
two ellipses. 


326. Proposed by Erich Michalup, Caracas, Venezuela. 


Sum the series 


SOLUTIONS 
Late Solutions 
293, 296. C. W. Trigg, Los Angeles City College. 
PARALLELEPIPEDS WITH INTEGER EDGES 
289. [Nov. 1956, May 1957] Proposed by J. M. Howell, Los Angeles City 
College. 


The diagonal of a rectangular parallelepiped can be any odd integer 
other than one or five and the edges integers relatively prime to the 
diagonal. 
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Comment by D. A. breault, Waltham, ass. 


The restated theorem, as was its predecessor, is false. \\e are 
considering solutions of 


x2 + y? + 22 = p2 (1) 
such that k is odd, not equal to | or 5, and 
(%,k) = (y,k) = (2,k) = 1. (2) 


Now for k = 15, the only positive integer solutions of (1) are 
(14, 5, 2), (11, 10, 2), and(10, 10, 5), no one of which satisfies (2). 


The contra-example disproves the proposition. 


AN INEQUALITY 
299. (March, 1957] Proposed byLeo Moser, | niversity of Alberta, Canada. 


te? prove that 


n!} < 


a'bict... atbec®,, 


Solution by Chih-yi hang, University of \innesota. With a minor 
change in notation, we write 


Each term of the multinomial expansion of the left side 1s positive 
and less than the sum of all the terms. Hence tlie particular term 


tk 
ns: 


The stated inequality follows immediately. The equality holds only 
if a, =e for positive integer values of n, 


Also solved by Leon Bankoff, Los Angeles, Calif.; E. S. Keeping, 
University of Alberta, Canada; G. B. Parrish, Durham, N. C.; D. A. 
Steinberg, University of California Kadiation Laboratory, Livermore, 


Calif.; and the Proposer. 


A VECTOR SUM 


300- (March, 1957] Proposed byH. wy. Gandhi, Lingraj College, be lgaum, 
India. 
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Three forces P, Q and k, taken in order, act along the sides of 
an equilateral triangle. These forces vary to become respectively 
L, M and N but their resultant remains the same in magnitude and 
direction. Prove that (#-N)/2 Q+M-L. 


Comment byC. W. Triggz, Los Angeles City College. Without restriction 
to an equilateral triangle, consider vectors P, Q, Kk acting along sides 
a, 6b, c, respectively, of a triangle. ThenP +O +R=L +& or 
Soh ~ N= 200 ~L) analy if 30 + = - 
But these vectors have different directions, so each must be nul! 
whence L = P/3 and M > -3Q0, whereupon N cannot act parallel to R in 
order to obtain the unchanged resultant. 

If the resultant is to remain unchanged it is necessary that the 
magnitudes of the forces be such that 


(L - P)/a = (M - Q)/b = (N - R)/c. 


Ek. S. Keeping, University of Alberta, Canada also concluded that 
the stated equality does not hold. 


A TRIGONOMETRIC IDENTITY 


301. (March 1957] Proposed by Norman Anning, Alhambra, California. 


Prove that 
Tan 5A =tan A tan (A + 36°) tan(A + 72°) tan (A + 108°) tan (A + 144°). 


I. Solution by Harley Flanders, University of California, berkeley. 


Let a = e27'/5, then 
(x5 - y5)/(x - y) = (ax - y) (a2x - y) (ax - y) (atx - y) and 
(= ¥) = (ex ¥) ¥)\ (09x + 7) + ¥)- 


In these two identities set x = e'A and y = e”'4. There result 
sin 5A/sinA = 16 sin (At7/5) sin (At2m/5) sin (At37/5) sin (At+477/5) 
cos 5A/cosA = 16 cos (At7/5) cos (A * 27/5) cos (At37/5) cos (At47/5). 

The stated relation is obtained by dividing the last of these 


identities into the previous one. 
I]. Solution by L. A. Ringenberg, Kastern Jllinois State College. 


ond 
It is well-known that tan*36° = 5 - and tan = 5 + NS. 

Now if one of the angles A, A + 36°, A+ 72°, At 108°, A + 144°, 
is an odd integral multiple of 7/2, then 5A is also an odd integral 
multiple of 7/2, and conversely. In this case each member of the 
given equation is undefined. 
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Suppose, then, that A is any angle such that none of the angles. 
A, A + 36°, A + 72°, A + 108°, A + 144°, is an odd integral multiple 
of 7/2. Then 


tan A tan(A + 36°) tan(A + 72°) tan(A + 108°) tan(A + 144°) 


= tan A tan(A + 36°) tan(A - 36°) tan(A + 72°) tan(A - 72°) 


tan _tan?A - - tan*72°_ 

1 1 - tan2A tan236° =] - tan272° 
_tanA - (5 - 245) - (5 + NS) 
1 - tan2A(5- 245) 1- tan’(5 + N5) 


_ tan tan4A - 0tan2A4 + 
Stan44 - 10tan2A + ] 


Using the formulas for tan(atf) and tan 2a we reduce 
tan 5A to the same expression. 4 


Q tan tanA(4 - 4tan2A) 
- tan224 tantA - 6tan2A + ] 


tan 4A = 


4tan2A 
__tan A + tan 4A tan 6tan’A 
tan 5A ~ tan A tan 4A Atan*A - 4tan4A4 

tan4A - 6tan2A + ] 


tantA - l0tan’A + 5 
Stan¥4 - +] 


Also solved by Leon Bankoff, Los Angeles; G. b. Parrish, Office 
of Ordnance Research, Durham, N. C.; C. W. Trigg, Los Angeles City 
College; Chih-yi Wang, University of Minnesota; and the Proposer. 


OVALS DEVELOPED FROM CURVES ON CYLINDER 


30z- [March 1957] Proposed by C. S. Ogilvy, Hamilton College, New York. 


If a 360° arc is drawn with a pair of compasses on the surface 
of a cylinder and the surface is then -unrolled« and laid flat, the 
resulting plane curve is an oval. 


(1) Characterize the oval if the cylinder is a right circular 


cylinder. 
(2) Characterize the cylinder if the oval is an ellipse. 


Solution by the Proposer. (1) The curve to be developed is the 
intersection of a sphere with a right circular cylinder, the center 


| 


116 MATHEMATICS MAGAZINE (Nov. -Dec. 


of the sphere lying on the cylinder’s surface. When unrolled the 
curve has the transcendental equation: 


x2 + d2 sin2(y/d) = r2 


where r is the radius of the sphere and d is the diameter of the 


cylinder. 

(2) Impossible. The cylinder must bea plane or the intersection 
of two planes which meet at the center of the compass arc, and of 
course in that case the ellipse is a circle. 


COMBINATORIAL SUMMATIONS 
303- [March 1957] Proposed by D. A. Steinberg, Livermore, Calif. 
Let gq, n be any positive, non-zero integers. Prove 
q(k -n) tk qtk 


qtk k 


- n2) + k(Qk - ‘lee 


. 
gtk k 


Solution by Chih-yi Wang, University of Minnesota. Mathematical 
induction will be used for the proofs. 
(1) is true for n = ]. Let us assume it is true for n = m. Then 


we have 


g(k-m-]) +k (2 + 
qt+k k 


wi k-m) +k qtk 
tk 


which is the required form for n =m + ]. 


€ 

k k=i gtk k qtmt+y] 

=m- + =m- +7" 

qtm q q +m 

= = + + = 
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(2) 1s true for n = ]. Let us assume it is true for n = m, Then 
we have, by aid of the argument in (]). 
q(k?-m?) + k(2k-1) gtk 
1 q tk k 


nti q{k2-(m+1)2] th(2k-1) 


= 


k=1 qtk k 


+ >. 


qtm 


qtm q 
) = m2-(Qm +1) - ( ) 


(m*1) (2m*+1) qtmt] 
+ 
=m2+ In + 1 = (m+ 1)? 
which 1s the required form for n =m + }, 


Also solved by Waleed A. Al-Salam, Duke University, N. C. 


CIRCLES CONNECTED WITH TRIANGLE 
304. [March 1957] Proposed by Huseyin Demir, Kandill, bolgesi, Turkey. 


Let AbC be a triangle, AB 7 AC, inscribed in a circle (O), and let 
K be the point where the exterior angle bisector of A meets (Q). A 
variablé circle with center at K meets Af, AC at E and F respectively, 
such that 4 is not an interior point of KkF. Find the limiting position 
m of the common point \ of kK, BC as EF approaches BC. 

Solution by the Proposer. Let E', F' be the points where (K) meets 
Ab, AC other than £,+#. Let M' be the common point of BC with E'F'. 

Applying the Menelaus theorem to ABC, considering EFM, E'F'M' as 
transversals, we have 


Multiplying these equalities member to member and observing that 
kA = F'A, E'A FA we get 


Since in the last ratio the numerator and denominator are the 
powers of 6, C with respect to the circle (K), and since these powers 


e qtk 
WB FC EAL ,, EYA 
FA’ EB wera’ 
OEB.E'B 
wc OFC.F'C 
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are equal (K is equidistant from 6 and C) WB:MU = M'C:M"B follows. 
Hence the points M and 4! are symmetric points on bC. The limiting 
position m of M will also be symmetric point of a!, the limiting 
position of M', [t is easy to see that m! is the foot of KA, the 
exterior angle bisector of A. Hence the construction of M follows 


immediately. 


LIMIT OF A SUM 


305- [March 1957] Proposed by y. S. Klamkin, Polytechnic Institute of 
brooklyn. 


Find gn 
im: fe * 2 
Solution by L. Cariitz, Luke University, N. C. Logarithmic 
differentiation of the familiar identity. 


For x = — this becomes 
z 


0 


and therefore 


Also solved by Chih-yi Wang, University of Minnesota; and the 


proposer. 


QUICKIES 


From time to time this department will publish problems which may be 
solved by laborious methods, but which with the proper insight may 


| 
0 
gives 
0 xé l-x 
2? | 
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be disposed of with dispatch. Readers are urged te submit their 
favorite problems of this type, together with the «Legant solution 
and the source, if known. 


Q@ 32. [March 1951] Each edge of a cube is a one-ohm resistor. What 
is the resistance between two diagonally opposite vertices of the 


cube? 


[Alternate solution by C. W#. Trigg.) Consider planes passed through 
the extremities of each group of three resistors issuing from two 
Opposite vertices. These planes divide the resistors into three 
series-connected sets of "parallel« resistors, consisting of 3,6, and 3 
resistors, respectively. Hence the total resistance is 1/3 + 1/6 + 1/3 
or 5/6 or 0.833 ohms. 

+1 n+ 


@ 211. Show that Os 2 ? is divisible by 7. 
(Submitted by yw. S. Klamkin) 


@ 212. When are the hour, minute, and second hands of a clock 
together? 

(Submitted by Barney Bissinger) 
Q 213 Prove that every skew symmetric determinant of odd order has 
value zero: 


(Submitted by M. S. Klamkin) 


Q 214. The three joins of the pairs of points of trisection of the 
sides of a triangle most remote from the vertices are concurrent. 


(Submitted by J. M. Howell) 


Q 215. Determine the minimum value of the sum of the squares of the 
 perpendiculars drawn from a point in the plane of a triangle to its 
three sides. 


‘(Submitted by S. Klamkin) 
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THE TREE OF MATHEMATICS 


The unusual nature of this book makes desirable some description in 
addition to the general information in the September-October issue of 
this magazine. 

The first four chapters, presupposing only arithmetic, introduce in a 
realistic style; algebra, plane geometry and trigonometry. The fifth chapter 
extends the natural number system (the integers) to include irrational 
numbers, thus completing preparation for analytic geometry and the calcu- 
lus. Each of these seven chapters concludes with a short list of exercises 
designed to illustrate and broaden the concepts which they have discussed. 

The succeeding twenty chapters are, for the most part, self-contained 
except that they presuppose the first seven, or the same subjects studied 
elsewhere. 

Puzzlement, formal manipulation and routine drill have been kept at a 
minimum. These matters play important psychological roles in conventional 
text books but would be out of place in this book whose purpose is to 
present the principles of mathematics as clearly and briefly as possible. 

Teachers who use this book for class work will probably want to supply 
exercises suited to their classes. 
suited to preparatory and general survey courses.) 

The kind of people who will use it for home study will, I believe, fix 
its discussions in mind by rereading and/or dreaming about them. 

Those who wish tostudy further into any of these topics can find detail- 
ed treatments in conventional textbooks on the classic subjects, while the 
chapters on newer subjects, e.g. ‘‘The Mathematical Theory of Games’’ 
and ‘‘Dynamic Programming,’’ contain extensive bibliographies. 

Finally, readers who do not care to concern themselves with technical 
terminology can get a good idea of the field of mathematics by reading 
just the expository introductions to the various chapters. 

The TREE OF MATHEMATICS when finished contained 420 pages 
instead of 350 as previously announced; but it still sells for $6, or $5.50 
cash with order. 

Friends of the Mathematics Magazine should know that it will recieve 
one half of the royalties from the Tree of Mathematics. 
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